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Ferromagnetism in purely organic materials is predicted for polymer and oligomers of dimethyl- 
aminopyrrole. First-principles calculations with the spin- density functional theory show that 
oligomers larger than tetramer as well as a polymer should have spin polarizations when the density 
of doped holes exceeds 1/2 hole per five-membered ring. We suggest that the mechanism for the 
ferromagnetism is, unlike the flat-band ferromagnetism predicted previously for polyaminotriazole, 
related to the periodic Anderson model, so that the present system amounts to an organic realization 
of the periodic Anderson model. 

PACS numbers: 71.20.Rv, 71.10.Fd, 75.50.Xx 



I. INTRODUCTION 

Realizing organic ferromagnets consisting entirely 
of non-magnetic elements is a challenging target in 
condensed-matter physics and chemistry. To stabilize 
spin polarization in a p-electron organic material, inter- 
atomic electron-electron interactions such as an intra- 
molecular coupling or an intra-polymer coupling can play 
an important role besides the intra-atomic direct and 
Hund's couplings. In usual organic materials, however, 
such interactions is usually not enough to guarantee spin 
polarizations. 

In an entirely different approach to overcome this dif- 
ficulty, the present authors have theoretically proposed 
to utilize the "flat-band ferromagnetism" , originally pro- 
posed by Mielke and by Tasaki^ to design an organic 
ferromagnet. This mechanism for ferromagnetism is dis- 
tinct from the usual band magnetism in that this is out- 
side the conventional Stoner picture, but instead arises 
from an electron correlation effect exerted on the flat (dis- 
persionless) band where the flat band is an outcome of a 
quantum mechanical interference rather than a narrow- 
band limit. However, a flat band has to satisfy what 
is called the local connectivity condition for the ferro- 
magnetism to be realized. This is a special condition 
in which the adjacent "Wannier" orbits have to overlap 
with each other no matter how they are combined to min- 
imize their size. Intuitively, the Bloch states on the flat 
band are something like non-bonding molecular orbitals 
in the language of organic chemistry, and the "general- 
ized Hund's coupling"— i.e., Hund's coupling in /c-space 
acts to realize the spin alignment when the condition is 
met. Since the flat-band model has eluded an unambigu- 
ous realization in real materials so far, we started from 
a materials-design viewpoint that it should be easier to 
satisfy the connectivity condition if we consider chains 
of odd-membered rings. So polymers of five-membered 
rings are candidates, but chains of odd-membered rings 



are by no means a sufficient condition for the flat band, so 
that we have examined a vast number of compounds. As 
a candidate for the realization of the magnetism, we have 
proposed polymers of five-membered aminotriazole^^^^^'^. 
Following the suggestion, Nishihara and coworkers have 
actually started fabricating the material.'' 

Given this background, the present paper has two 
strong motives that have led us to a new avenue. Exper- 
imentally, it has turned out that polymerizing monomers 
of aminotriazole, which is commercially available, has 
turned out to be difficult, presumably due to a quantum 
chemical reason. In the course of our attempts, how- 
ever, we have found that monomers of dimethylaminopyr- 
role can be polymerized more easily, where triazole part, 
presumably the difficult ingredient, is absent while the 
amino group is replaced with dimethylamino group. ^ 
Thus we have started with obtaining the electronic struc- 
ture of polymer and oligomers of dimethylaminopyrrole 
from first principles to examine the possibility of ferro- 
magnetism (for the polymer) and high-spin states (for the 
oligomers). We have employed the generalized gradient 
approximation (GGA)^. 

New findings here are two-fold: we shall first show for 
the polymer that the materials will be indeed spin po- 
larized if holes are sufficiently doped (i.e., more than 1/2 
hole per five-membered ring). Oligomers are also shown 
to have high-spin states if holes are doped enough. For Ti- 
mers of dimethylaminopyrrole, we conclude that n/2 + 2 
or n/2 + 3 holes are necessary to be ferromagnetic de- 
pending on whether n/2 is even or odd, respectively. 

More importantly, from the viewpoint of mechanism of 
ferromagnetism, the second finding is that the band on 
which the spin polarization resides does not satisfy the lo- 
cal connectivity condition although it is quite dispersion- 
less, so that the mechanism of the ferromagnetism should 
indeed be outside the flat-band ferromagnetism. We shall 
then discuss that the magnetism can be explained by the 
mechanism proposed by Batista et al.^ based on the pe- 
riodic Anderson model. In this theory, coexistence of a 
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FIG. 1: A top view(a) and side view(b) of the optimized 
atomic structure of polydimethylaminopyrrole in a unit cell, 
(c) Band structure of undoped polydimethylaminopyrrole. 
Thick solid (dashed) lines represent cr-bands (vr-bands). 

narrow and correlated band hybridized with a dispersive 
and uncorrelated band is a key to ferromagnetism. Be- 
cause that mechanism is a kind of generalization of the 
atomic Hund's rule to the solid, it is natural that spin- 
polarized states are obtained by DFT calculations. If this 
picture is correct, the present prediction for ferromag- 
netism in dimethylaminopyrrole polymer and oligomers 
amounts to a novel realization of the periodic Anderson 
model in purely organic materials. 

II. DIMETHYLAMINOPYRROLE POLYMER 

Let us start with the polymer of dimethylaminopyr- 
role. Figure [Ha,b) depicts the atomic structure of poly- 
dimethylaminopyrrole in a unit cell. We have studied 



both the spin-unpolarized cases based on the density- 
functional theory (GGA-DFT) and the spin-polarized 
ones based on the spin-density functional theory (GGA- 
SDFT). We have used the plane- wave based ultrasoft 
pseudopotentials^^'^^ with the energy cutoff at 20.25 Ry- 
dberg. The convergence criterion of the geometry opti- 
mization was that all of the forces acting on each atom 
were within 1 x 10~^ H/a.u. We have first optimized the 
structure. While we have started, in our materials de- 
sign, from a planar configuration of rings, the structure- 
optimized polydimethylaminopyrrole has dimethylamino 
bases somewhat twisted out of the plane. This makes 
the TT-orbital of the nitrogen atom in the amino base ly- 
ing in the molecular plane, that is, the angle between 
the TT-orbitals of the nitrogen atom in the amino base 
and the neighboring nitrogen atom in the five membered 
ring becomes nearly 90 degrees, leading to a very small 
TT-bonding between them. 

The band structure of undoped polydimethylaminopy- 
rrole is shown in Fig. HJc). We can see that there exist 
flat bands (labeled as (B)-(E) and (C)-(F) in Fig.[I{c)), 
but they lie well below the Fermi energy, unlike the case 
of polyaminotriazole^i^i^ where the flat band sits around 
Ef- In the band structure, a band folding of the quasi- 
one-dimensional band is due to the alternating directions 
of the pyrrole molecules. 

Figure [2](A-H) shows the wave functions in the high- 
est four valence (occupied) bands at F and X points. 
The wave functions of the flat bands (B, C, E and F in 
Fig. [IJc)) are mixtures of the 7r-orbitals of the nitrogen 
atoms and a-orbitals of the carbon atoms, while the other 
wave functions consist almost entirely of 7r-orbitals. We 
can thus call the flat bands (B)-(E) and (C)-(F) (thick 
solid lines in Fig. W^c)) as <j-bands, and the top valence 
band (A)-(G) (uppermost thick dashed line in Fig. [T]) as 
a TT-band. 

While in polyaminotriazole^^^ the 7r-band is flat, for 
which the "flat-band ferromagnetism" is expected to be 
realized, this is no longer the case with polydimethy- 
laminopyrrole, since the 7r-band is considerably disper- 
sive. On the other hand the cr-bands are flat, and one 
might expect that the fiat-band ferromagnetism may oc- 
cur. This does not apply either, since the bands do 
not satisfy the local connectivity condition^ mentioned 
above, but they are fiat simply because the transfer be- 
tween five- membered rings is small, as seen in Figs.[2fB), 
(C), (E) and (F) where a-bonds are absent between those 
two carbon atoms connecting the adjacent five- membered 
rings. Although small 7r-like orbitals can be seen there, 
they are negligibly small compared to those of 7r-bands 
in Figs.El^G) and (H). 

Does this imply that we have no ferromagnetism? If 
we look at the band structures of the hole-doped poly- 
dimethylaminopyrrole calculated with GGA-SDFT for 
various values of the hole concentration in Fig. [3l while 
there is no spin polarization for the number of holes per 
unit cell nh = 0.5 or 1.0, spin polarization does appear 
when nh exceeds 1.0. The differences between the num- 
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(a) nh = 0.5 



(b) nh=1.0 



FIG. 2: Wave functions in the highest four valence bands of 
undoped polydimethylaminopyrrole at F (A-D) and X (E-H) 
as depicted by isosurfaces with different colors representing 
the signs of the wave function. See Fig. [T]^c) for labels of the 
energy levels. We have chosen the wave function to be real at 



bers of up and down spins per unit cell are 0.17, 0.45 and 
0.97 for nh = 1.25, 1.5, 2.0, respectively. This shows that 
more than half- filled holes (nh > 1.0) make the system 
ferromagnetic. This is the first key result. 



III. ORIGIN OF THE MAGNETISM: PERIODIC 
ANDERSON MODEL 

So the natural question is what is the mechanism of 
the spin polarization observed in the present calculation, 
which should be other than the flat-band mechanism. 
We can evoke ferromagnetism in the periodic Anderson 
model as a possible scenario. Indeed, here we have nar- 
row (T bands embedded in wide tt bands which hybridize 
with each other. Theoretically, it has been well known 
that the periodic Anderson model can have a ferromag- 
netic ground state in a variety of conditions'^. One fa- 
mous example is the situation when the model can be 
mapped to the Kondo lattice model^^, i.e., when there 
is one electron in each localized orbital. On the other 
hand, Batista et al.^ recently proposed another mecha- 
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FIG. 3: Spin-resolved band structures of the hole-doped poly- 
dimethylaminopyrrole for the concentration of holes per unit 
cell (a) rih = 0.5, (b) rih = 1.0, (c) uh = 1.5, (d) rih = 2.0. 



nism which can work when the particle density of the 
localized state is not close to unity. 

The basic idea of Batista et al^ may be recapitulated 
as follows. Let us consider the periodic Anderson model 
whose Hamiltonian is defined as 
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where the conduction electrons (created by c}, ^) hy- 
bridize with localized electrons (created by /^^) with r 
representing the site and a the spin. Here tc^tf are the 
respective hopping energies, Sf (assumed to be \ef\ < tc) 
the energy level of the localized orbital, V the hybridiza- 
tion, and U is the repulsive interaction within the local- 
ized orbital. We can diagonalize the one-body part of the 
two-band Hamiltonian Hq to obtain the band represen- 
tation as 



(4) 



where > (see Fig El . 



Let us introduce two kinds of Wannier orbitals, one 
dominantly localized, 

K> K< 

k k 

and another dominantly conductive, 

K< K< 

k k 

where we have decomposed the /c-space into {K^) for 
which the localized (conduction) electrons are the domi- 
nant in the Pk,a band (FigE]). 

With the new Wannier basis the one-body part of the 
Hamiltonian now reads 

where r's are new transfer energies with no off-diagonal 
terms like il)l^^(j)r' ^a- Batista et al.^ have shown that, 
while the interaction term in the Hamiltonian becomes 
complicated in terms of ip and 0, the dominant term is 
Uef^nf^nf, while the other terms can be neglected when 
the hybridization V is not too strong. When the system is 
in the mixed valence regime, namely, when the Fermi en- 
ergy Ef is above Sf in the present hole-doped case, holes 
first doubly-occupy the uncorrelated i/jk states which are 
above £/, and then go into (j)k states. There (j)kS tend 
to have parallel spins, because the cost of the kinetic en- 
ergy due to the spin polarization (~ tf) is small, while 
the ferromagnetic state does not feel /7efr because of the 
Pauli's exclusion principle. 

When all the ^^'s are singly occupied (by holes in the 
present case), we can construct localized 0r's by using 
occupied states only. Because there is no interaction 
between the spins in this case, the ferromagnetic state 
and the paramagnetic state have the same energy. When 
0/c's are only partly occupied, by contrast, we cannot 
construct completely localized (j)r^s only from the occu- 
pied states. Carriers then start to feel Uef^ through the 
overlap of delocalized 0r's, which makes the ferromag- 
netic state energetically favorable. Therefore, the sys- 
tem should choose a ferromagnetic ground state when (j)k 
states are less than half-filled for [/gff sufficiently greater 
than t f . This has been demonstrated numerically for the 
one-dimensional periodic Anderson model*^ 

In this mechanism of ferromagnetism, a hybridization 
between narrow and dispersive bands plays an important 
role. In polydimethylaminopyrrole, it is clear from the 
band structures in Fig. [U^c) and Fig|3] (which are sum- 
marized in Fig|5]) the a-band and 7r-band are hybridized 
each other. Although the two bands may seem intersect- 
ing, they are actually anticrossing with each other. So if 
this scenario is correct, the organic material considered 
in the present work amounts to the first proposal for an 
organic periodic Anderson system. 
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FIG. 4: Band structure of polydimethylaminopyrrole in a 
tight-binding model. The transfer energies and on-site en- 
ergies are fitted to reproduce Fig. [T]^c) . Thick solid (dashed) 
lines correspond to E'^ (E^) in Eq. Q. (A-D) on the right 
depict wave functions for the highest four occupied states at 
r, where the size of each circle represents the amplitude while 
the different colors its sign. 



IV. TIGHT-BINDING ANALYSIS 

In order to grasp what is essentially happening in the 
ferromagnetism of this material, we introduce a simple 
tight-binding model. The model consists of five sites rep- 
resenting the TT-orbitals in the five-membered ring and an 
additional "amino N" site. The amino N site represents 
the molecular orbital that consists of the 7r-orbital of the 
N atom in the amino base, the a (sp^) orbitals of two 
methyl bases and a-orbitals of the five-membered ring. 
The transfer and on-site energies are fitted to reproduce 
Fig. [T](c). Figure m shows the band structure and wave 
functions in the tight-binding model. In the tight-binding 
model the primitive cell contains only one five-membered 
ring, but we have folded the bands at X (/c = 7r/2a) to 
facilitate comparison with Fig. 1(c). We can see that the 
tight-binding fit successfully reproduces all the relevant 
(i.e., the thick lines in the DFT result in Fig. 1(c)) tt- 
bands and a-band. So the six sites considered in the 
tight-binding model have turned out to be sufficient for 
describing the relevant bands (consisting of the mainly 
TT bands of the ring, plus another one whose character is 
not the TT orbits of the ring). In other words, the bonding 
between the ring and amino N can be expressed only by 
the transfer between the 7r-orbital of the N atom at the 
apex of the ring (apex N) and that at the amino N. Be- 
cause the molecular orbital that forms a-band is included 
in the amino N site, the a-band of this model has a large 
amplitude at the amino N. In this sense, "a-band" (which 
was meant to mean those other than the tt orbits of the 
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FIG. 5: Band structure of polydimethylaminopyrrole in the 
tight-binding model, unfolded onto an extended Brillouin 
zone. Solid (dashed) lines correspond to E'^ (E^) in Eq. (|4]), 
while thick (thin) lines ('0) states in Eq. ^ a (tt) bands 
in the present system). Labels (A-D) refer to the band struc- 
ture in Fig. U Arrows schematically indicate the spins of 
accommodated holes. 



ring) should rather be called "dimethylamino band" , or, 
put more generally, the functional group band. 

The simple model shows that the a-band does not sat- 
isfy local connectivity condition. The flatness of the band 
simply comes from localized nature of the wave function. 
While the tight-binding model here is the same as that 
for polyaminotriazole^ except for the values of parame- 
ters, a large difference is that the transfer between the 
amino N and apex N is very small here compared to that 
in polyaminotriazole, which is precisely because the tt- 
orbitals of the N atoms are twisted almost 90 degrees 
with each other here. The twist also produces hybridiza- 
tion with cr-orbitals of the five-membered ring. The small 
transfer between the amino N and apex N breaks the 
subtle balance of parameters necessary to satisfy local 
connectivity condition for the flat-band ferromagnetism, 
but simultaneously creates narrow and dispersive bands, 
mutually hybridized, which is a prerequisite for the fer- 
romagnetism based on the periodic Anderson model. 

In order to make the identification of the present band 
structure with the periodic Anderson model clearer in 
terms of (j) and i/j states in Eq. (j6|), we show in Fig. [5] 
the band structure (Fig. [4j) unfolded into an extended 
Brillouin zone. The upper band (solid line) consisting 
of a part of the top valence 7r-band and a part of the 
flat cr-band correspond to in Eq. (j3j) , while the lower 
band (dashed line) consisting of the other parts of a and 
TT-bands correspond to . In the figure, thick lines and 
thin lines correspond to (j) and ?/^-states in Eq. ([6]), re- 
spectively. Simply put, (j) represents the a-band and 
the TT-band in the present organic system. 

In terms of the tight-binding model, we can explain 
an interesting feature in the GGA-SDFT band structure 
(Fig. [3|) where the energy of the a-band shifts relative to 
the other bands with the doping. For instance, as we in- 
crease the number of holes from nh = (Fig. [31(a)) to 1.0 
(Fig. [3](b)), the a-band goes up. When nh is increased 
further (Figs. [3l^c) and [31(d)) , the a-band of majority 
spins starts to go down, while that of minority spins does 
not. As a result, the lower limit of doping concentra- 
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FIG. 6: Atomic structure(a) and energy levels of undoped(b) 
and four- hole-doped (c) tetramer of dimethylaminopyrrole. 
SoUd (dashed) lines represent occupied (unoccupied) states. 
Wave functions of the states labeled as (A-L) are shown in 
Fig. [3 below. 



tion necessary for the ferromagnetism becomes just half- 
filled (nh = 1.0, i.e., 1/2 hole per five membered ring). 
Such a shift of the <j-band in the GGA-SDFT calculation 
can be reproduced by using the tight-binding model, if a 
Coulomb interaction of 1.4 eV between the amino N site 
and five-membered ring is assumed and a mean-field ap- 
proximation is applied. Note that tt orbital of the amino 
N is hybridized with the a orbitals of the five-membered 
ring, so that the interaction effectively contains an intra- 
site interaction and can be large. 



V. OLIGOMERS 

We move on to the oligomers of dimethylaminopyr- 
role. This is practically important because synthesis of 
the oligomers of dimethylaminopyrrole is possible^ while 
the polymer may be difficult to fabricate. 

As a typical case we show the result for a tetramer. 
Figure [6] shows the atomic structure, the energy levels of 
undoped and four- hole-doped tetramers. For the hole- 
doped tetramer, we can see in Fig. [6](c) that four con- 
secutive minority-spin states (H-K) take relatively higher 
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FIG. 7: Wave functions of a four- hole-doped tetramer of 
dimethylaminopyrrole for the energy levels labeled in Fig.fGlJc) 
as depicted by isosurfaces with different colors representing 
the signs of the wave function. 



energy levels than corresponding majority-spin states (C- 
F), while other levels are almost the same between the 
majority and minority spins. As a consequence the higher 
two states out of the four become unoccupied in the 
minority-spin channel. 

The four states are <j-states, counterparts of the <j- 
band in the polymer, which can be confirmed if we look 
at the wave functions of the four states in Figs. [7l(C-F) 
for majority spin and in Figs. [Tl^H-K) for minority spin. 
All these wave functions are seen to consist of nitrogen 
TT-orbitals and carbon a-orbitals, while the other wave 
functions in Figs. [7l(A,B,G,L) consist only of 7r-orbitals. 

We can relate the arrangement of energy levels of the 
oligomer with the band structure of the polymer. Gener- 
ally, for an oligomer consisting of n monomers (n-mer), 
the number of electronic states constructed from one kind 
of molecular orbital of the monomer is n. Obviously, a 
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FIG. 8: Correspondence between the band structure of the 
polymer (curves) and the energy levels of oligomers (circles) 
for dimethylaminopyrrole for the four-hole-doped tetramer 
(a) and the six- hole- doped hexamer (b). Single, double, and 
triple circles indicate one, two and three energy levels in the 
oligomer, respectively. 



state in which all monomers have the same phase is one 
of the n states, and its wave function is similar to that 
of the polymer at F-point except for the truncation of 
the chain. We call it F-state. Another linearly inde- 
pendent state can be constructed by introducing a node 
into the wave function of the F-state making the signs 
of the molecular orbitals at the opposite ends reversed. 
This state is similar to that of the polymer at /c = 7r/na 
(a: the length of the monomer) in that the phase ro- 
tates by TT per n monomers. In the same way, m-th state 
of the oligomer corresponding to the polymer's state at 
k = (m — l)7r/na can be constructed by introducing m — 1 
nodes (a phase rotation of (m — l)7r) in the F-state. Be- 
cause of the similarity of the wave functions, the energy 
level of the m-th state in n-mer can be inferred according 
to the band energy of the polymer at /c = (m — l)7r/na. 

An easy way to do such an estimation is to take the in- 
tersecting points of the band with vertical lines dividing 
the zone in k-space {k = to ir/a) into n. For example, 
the energy levels of the four-hole-doped tetramer and six- 
hole-doped hexamer around the LUMO-HOMO can be 
estimated as in Figs. [51^ a) and[Hb), respectively. Note 
again that the whole band is folded at X (/c = 7r/2a) be- 
cause of two monomers in the unit cell. The band struc- 
ture of Fig. [3l^d) is used here because the hole density 
of these examples is nh = 2.0. We take n to be even in 
the following for simplicity. As shown in the following, 
the order of energy levels estimated in this way agrees 
well with the result for the oligomers. One reason why 
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TABLE I: Difference between the number of up and down spins are shown for various ohgomers (rows) with various numbers 
of holes (columns). The numbers in the parentheses indicate the energies of the spin polarized states (GGA-SDFT) measured 
from those for unpolarized states (GGA-DFT). 
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6-mer 


0.00 


0.00 
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8-mer 


0.00 


0.00 


1.84 (-0.07) 


3.22 (-0.33) 






10-mer 


0.00 


0.00 


0.32 (-0.03) 


2.09 (-0.07) 


4.07 (-0.46) 




1 2-mer 


0.00 


0.00 


0.00 


1.84 (-0.01) 


3.06 (-0.16) 


4.97 (-0.57) 



such an estimation is valid is that the finite-size effect is 
not very large in the oligomers of dimethylaminopyrrole, 
probably due to the simple chain structure and not-too- 
large transfer between monomers. 

With such a kind of level scheme, we can estimate how 
many holes should be doped to realize spin polarization in 
a given oligomer. For that purpose, whether the vertical 
line in Fig. [8] exist at k =X/2 is important, since the tt- 
band and the flat <j-band of minority spin cross at this 
point. When n/2 is odd as in Fig. [Hl^b), no vertical line 
exists at X/2 and thus n/2 + 1 holes first enter the levels 
coming from the 7r-band whose energies are the same for 
both spins. Since two more holes are necessary to occupy 
the levels associated with the flat a-band, n/2 + 3 holes 
are necessary in total for the state to be spin-polarized. 
The total spin (i.e., the difference, r^, of the number of 
up and down spins) can be estimated to be A/" — n/2 — 1, 
where N is the number of doped holes in an n-mer. When 
n/2 is even, on the other hand, one of the vertical lines 
passes through X/2 as in Fig. [81(a). At k = X/2, the 
four energy levels, two from the 7r-bands for both spins 
and two from the <j-bands for the minority spin, almost 
coincide when holes are doped. Depending on which level 
is the highest, n/2 + 2orn/2 + 4 holes are necessary to 
have spin-polarizations in this case, with rj estimated to 
be — n/2 or X — n/2 — 2, respectively. 

TableUshows the result of GGA-SDFT calculations for 
a series of oligomers of dimethylaminopyrrole with vari- 
ous numbers of doped holes. The number of holes nec- 
essary for spin polarization almost fully agree with the 
above estimate in terms of the polymer's band structure. 
The only exception in this table is the existence of a spin 
polarization in the six-hole doped decamer, but the spin 
polarization of 0.32 there is rather small. For the n-mers 
with odd n/2, i.e., hexamer and decamer, the values of 
the spin polarization agree accurately with the above es- 
timate, which should reflect the fact that no energy levels 
associated with the 7r-band exist near the Fermi energy 
when spin polarization emerges, as shown in Fig. [Sjb). 
For the oligomers with even n/2, i.e., octamer and dode- 
camer, the spin polarization takes the value somewhere 
between the two estimations, N — n/2 and N — n/2 — 2. 
In reality, 77 is expected to take either of those two values, 
depending on the order of almost degenerated energy lev- 



els of TT-band and a-band. The calculated values between 
the two estimations should come from fractional occupa- 
tions of these energy levels on the assumption of a finite 
temperature introduced to ensure the convergence of the 
electronic states. 



VI. CONCLUSION 

We have used a first-principles calculation to exam- 
ine magnetism in a polymer and oligomers of dimethy- 
laminopyrrole from a viewpoint of materials design for 
ferromagnetism without magnetic atoms. We have found 
that polydimethylaminopyrrole should be ferromagnetic 
if density of doped holes are more than one hole per unit 
cell {nh > 1.0). For n-mers of dimethylaminopyrrole, we 
conclude that n/2 + 2orn/2 + 3 holes are necessary to be 
spin-polarized depending on whether n/2 is even or odd, 
respectively. While polydimethylaminopyrrole belongs to 
the same category (chains of five-membered rings) as the 
previously proposed polyaminotriazole^^^ which has theo- 
retically shown to realize the fiat-band ferromagnetisn>i, 
we suggest that the present ferromagnetism is distinct 
from the fiat-band ferromagnetism, but has instead the 
mechanism proposed by Batista et al.^ based on the pe- 
riodic Anderson model. 

So the present work gives an interesting instance of ma- 
terial design where control of the parameters of electronic 
states even leads to a different mechanism for the many- 
body effects, namely organic periodic Anderson model 
and associated ferromagnetism. An important future 
problem is whether the SDFT calculations can capture 
the ferromagnetism d la Batista. Another future problem 
is how the present concept of the organic periodic An- 
derson model can be extended to wider materials beyond 
those considered here. As for the doping the material, we 
can consider various chemical dopants as discussed for in 
Ref. [i]. 
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